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ABSTRACT 

We report initial results from the far-infrared fine structure line observations of a sample of 44 local 
starbursts, Seyfert galaxies and infrared luminous galaxies obtained with the PACS spectrometer on 
board Herschel. We show that the ratio between the far-infrared luminosity and the molecular gas 
mass, Lfib./ M}i^, is a much better proxy for the relative brightness of the far-infrared lines than Lfir 
alone. Galaxies with high Lfir/A/hj ratios tend to have weaker fine structure lines relative to their 
far-infrared continuum than galaxies with Lfir/-M^H2 ^ 80 Lq Mq~^ . A deficit of the [C II] 158 /xm 
line relative to Lfir was previously found with the ISO satellite, but now we show for the first time 
that this is a general aspect of all far-infrared fine structure lines, regardless of their origin in the 
ionized or neutral phase of the interstellar medium. The Lfir/-^H2 value where these line deficits 
start to manifest is similar to the limit that separates between the two modes of star formation 
recently found in galaxies on the basis of studies of their gas-star formation relations. Our finding 
that the properties of the interstellar medium are also significantly different in these regimes provides 
independent support for the different star forming relations in normal disk galaxies and major merger 
systems. We use the spectral synthesis code Cloudy to model the emission of the lines. The expected 
increase of the ionization parameter with Lfir/Mhj can simultaneously explain the line deficits in 
the [C II], [Nil] and [O I] lines. 

Subject headings: galaxies: evolution — galaxies: ISM — galaxies: starburst — infrared: ISM 



1. INTRODUCTION 

During most of their life, galaxies form stars at a rela- 
tively low rate, regulated by t he total amount of molec- 
ular gas mass available (e.g.. iBigiel et al.ll2008D . This 
phase of semi-continuous star formation can be inter- 
rupted by short periods of violent star formation, nor- 
mally associated with the merging of massive galaxies. 
[ In mergers most of the molecular gas is compressed on a 
scale of a few hundred parsec to a few kpc around the nu- 
clei and, by a mechanism not yet well understood, forms 
[ stars with high efficiency (Sanders & Mirabel 1996). 
These episodes of intense star formation do not cover 
the whole merging process and might not be activated 
' in all types of galaxy interactions, since many optically 
classified mergers do not show at present a significant in- 
■ crease in their star formation rates (e.g.. iBushouse et al.l 
' [19881). 

The ratio between the far-infrared luminosity and 
the molecular gas mass, _Lfir/-^^H2i is expected to 
be proportional to the number of stars formed in 
the galaxy per unit molecular gas mass and time 
(Mh2 includes a 36% mass correction for helium, 
and Lfir is assumed to be star formation domi- 
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nated). Molecular line observations in the local Uni- 
verse have shown that on average galaxy mergers tend 
to have higher Lfir/-Wh2 ratios (> 100 Lq Mq~^) than 
non- or weakly-interacting star-forming galaxies (e.g., 
lYoung et al1ll986l: iSolomon fc Sagjll988f ). In addition, 
their far-infrared spectral energy distributions are con- 
siderably warmer (IRAS Seojim/ Sioo fj,m colors higher 
than 0.6: iMazzarella et al.lll99ll ). probably indicating a 
dust temperature increase due to the intense UV radi- 
ation fields associated with the high concentration of 
young stars. 

The effects of merger-induced versus quiescent star for- 
mation should also be reflected in the emission line prop- 
erties of the star-forming gas. For example, the intensi- 
ties of some fine structure lines produced in HII regions 
and photodissociation regions (PDRs) are strongly de- 
pendent on the value of the ionization parameter at the 
surface of the clouds (e.g., Abel et al. 2005), where the 
ionization parameter U is defined as the number of hy- 
drogen ionizing photons {hv > 13.6 eV) per hydrogen 
particle. As a first order approximation U should be pro- 
portional to Lfir,/-^H2 because the number of massive 
stars (and hydrogen ionizing photons) per unit of molec- 
ular gas mass increases with this ratio. High ionization 
parameter values have been proposed to explain the low 
[C II] 158 fim line to far- infrared continuum (FIR0) ratios 
observed with the ISO satellite in some luminous and 
ultr aluminous infrared galaxies (LIRGs a nd ULIRGs; 
e.g., iLuhman et all 120031: lAbel et al.l l2009f ). There are 

^ Wc use the definition of the far-infrared flux, FIR(42.5 /xm- 
122.5 /im), and the far-infrared luminosity, LpiR (40 /^m-500 /im) , 
given in table 1 of . Sanders fc Mirabel (1996 ). 
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Figure 1. The global [C II] 158 /^m line to FIR continuum ratio in galaxies with different optical activity classifications as a function of 
^FIR (a) and Lfir/Mhj (b). PACS observations are represented with squares. Open symbols indicate 3(t upper limits to the line flux. 
Grey areas connect points at one standard deviation below and above the mean (black line) for various Lpip and iFiR/AfH2 bins. 



however other possible explanations to this "[C II] line 
deficit" (saturatio n of the [C II] line, opt ical depth effects, 
soft UV fields; see lLuhman et aril2003 l for a detailed de- 
scription), that could not be fully excluded in the past 
due to the small number of infrared galaxies detected in 
far- infrared fine structure lines other than [C II] . 

Here we report initial results from the far-infrared fine 
structure line observatio ns of a sample of loc al galaxies 
obtained with Herschel (jPilbratt et al.ll2010( ). We find 
that galaxies with LpiR/Afna ^ 80 Lq Mq~^ tend to 
have lower line to FIR continuum ratios than galaxies 
with lower Lpm^/Mn^ values. These line deficits affect 
all the observed lines, independent of their origin in the 
ionized or neutral phase of the interstellar medium. 

2. OBSERVATIONS AND DATA REDUCTION 

The observations discussed in this Letter are part of the 
Herschel guaranteed time key program SHINING (P.I. 
E. Sturm), devoted to study the far- infrared properties 
of a sample of more than 100 galaxies, that includes lo- 
cal starbursts, Seyfert galaxies, low-metallicity systems 
and infrared luminous galaxies at low and high redshift. 
Here we present the fine structure line data for the first 
44 galaxies: 16 Seyferts (NGC 1068 and Gen A, among 
others), 5 HII galaxies (e.g., M 82, M 83, NGG 253), 21 
LIRGs and ULIRGs (e.g., NGG 4418, Arp 220, Mrk 231) 
and two high redshift galaxies (IRAS F10214+4724 and 
MIPS J142824.0+352619). 

The observations wer e obtained with the P AGS far- 
infrared spectrometer (jPoglitsch et al.l I2010D between 
October 2009 and December 2010. We observed the three 
PDR hues [G II] 158 ^Im, [O I] 145 ^J.m and [O I] 63 ^m, 
and the four HII lines [N II] 122 ^m, [O III] 88 ^m, 
[N III] 57 /im and [O III] 52 /im in most galaxies, with 
some exceptions due to data acquisition problems and 
scheduling constraints. The rest-frame wavelength of the 
[O III] 52 /im line is at the low sensitivity edge of the 
PAGS range and was observed only in the high-redshift 
galaxies (see lSturm et al.i r 20 10'). 

The data were reduced using the standard PAGS re- 
duction and calibration pipeline included in HIPE 5.0. 
However, for the final calibration we normalized the spec- 
tra to the telescope flux and recalibrated it with a ref- 
erence telescope spectrum obtained from dedicated Nep- 
tune observations. For unresolved sources, this method 



gives excellent agreement (differences <15%) between 
the IRAS 60 /im and 100 fim flux densities and those 
measured with PAGS when extrapolating from the con- 
tinuum close to the lines. 

Most of the targeted lines were detected. For the ex- 
tended galaxies, we integrated the line fluxes and FIR 
continuum within the 47"x47" PAGS field of view. The 
detailed spatial information for these sources will be dis- 
cussed in a future paper. 

3. RESULTS 

3.1. The [ClI] Line Deficit 

The well known [G II] line deficit in infrared luminous 
galaxies is shown in Figure [T^. Together with the PAGS 
data, we have plotted a compariso n sample coinposed 
of all galaxies observed with ISO (INegishi et al.l [20011: 
Malhotra et al.ll2001l:lLuhman et alll2003l: IBrauher et al.l 



2008t ). High redshift galaxies (z = 1.1-6.4) with 
[G I I] observations are also included (iMaiolino et alj 
'2009': 'Hailev-Dunsheat h et all 120101: ITyison et al.l 120101 
Wagg ct al. 2010; Stac ev et al.l I2010[ and references 
therein). The average [G II] /FIR ratio is approximately 
constant in galaxies with low infrared luminosities, but 
starts to decrease at Lfir ^ 10^^ Lq. At the same time, 
the scatter in the ratio increases by almost a factor of 
two: some high redshift galaxies do not show a [G II] 
deficit, even if their far-infrared luminosities are higher 
than 10^^ Lq, and the source with the lowest [G II] /FIR 
ratio is NGG 4418, a galaxy with Lfir — lO" Lq. 

The scatter at high Lfir can be reduced by a fac- 
tor ~1.5 if Lfir/Mhj is usecQ instead of Lfir (Fig- 
ure [TJa), presumably because this ratio is more closely 
related to the properties of the clouds (gas density and 
dust temperature) than Lfir alone (e.g.. lGao fc SolomonI 
120041 ). Galaxies with similar Lfir luminosities can 
have very different Lfir/-Mh2 ratios and the excitation 
conditions in their interstellar medium will also differ 

^ AfH2 was estimated using previously published CO obser- 
vations of the galaxies. We used a Xqq conversion factor 
that is a continuous function of the IRAS Sao fira/ Sioo iim color 
in a way that produces adequate values for Milky-way type 
galaxies {'^4MqL'~^, with L' in units of Kkms-^pc^) and 
ULIRGs (~O.8M0 L'-l). Some high redshift galaxies did not 
have enough photometric data points to estimate their rest-frame 
•S'so /im/S'ioo /im color; Xqq = 2 Mq L'~^ was used in those cases. 
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Figure 2. Global line to FIR continuum ratios in galaxies with different optical activity classifications as a function of Lpjpi/Muj ■ PACS 
observations are represented with squares. Open symbols indicate 3cr upper limits to the line flux. All the line to FIR ratios start to 
decrease at Lpijj/Mhj ~ 80 Lq Mq~^ (vertical lines). The horizontal lines indicate the approximate minimum values for the different line 
to FIR ratios observed in galaxies with Lpjr/Mhj < 80 Mq~^ 



(|Gonzalez- Alfonso et aLll2008D . With this new represen- 
tation, high redshift galaxies overlap with their low red- 
shift counterparts and NGC 4418 is shifted to the right of 
the figure due to its relatively low molecular gas c o ntent . 
A similar result was obtained by iMalhotra et"all ()2001f ) 
using the IRAS S'eo/Lim/'S'ioo/im color in a smaller sample 
of local galaxies observed with ISO. 

3.2. A General Far-infrared Line Deficit 

Figure [2] shows the line to FIR continuum ratio against 
Lfir/ for all the fine structure lines observed with 
PACS. Our main result is that we now find line deficits in 
both the FDR and the HII lines. All the line to FIR ra- 
tios start to decrease at Lfib^/M^^ ~ 80 Lq Mq~^. This 
result also does not seem to depend on the optical ac- 
tivity classification of the galaxies. We note, however, 
that the integrated emission of the [O I] and [O III] lines 
in some Seyfert galaxies might have a significant contri- 
bution from X-ray dominated region s close to the active 
galactic nucleus (AGN; e.g.. Malonev et al.lll996( ). 

The [O I] 63 fj,m line was detected in absorption at 
some velocities in NGC 4418, Arp 220 and IRAS 17208- 
0014. The [O I]63/im/FIR ratios used in Figure H cor- 
respond to the emission components only and should be 
regarded as lower limits. This will reduce the deficit ob- 
served for this line to some extent, b ut will not eliminate 
it, as other galaxies (e.g., Mrk 231: iFischer et al.ll2010( ) 
with no clear absorption features in their line profiles 
do still have lower [O I] 63 //m/FIR ratios than galaxies 
with Lfir/A/h2 < SOLqMq-'^. The [O I] 145 fim line is 
less likely to be affected by this kind of absorption and 
also shows a deficit at high LpmjMu_2- There is no sign 
of absorption by foreground material in the other fine 
structure lines at the current PACS velocity resolution 
(~100-250kms-i). 



In summary, our observations show that the line to 
FIR continuum ratio drops by a factor of 3 to 10 above 
Lfir/-^H2 80LqMq~^ in all the far-infrared fine 
structure lines we have observed. This drop seems to 
be a universal feature of galaxies at different redshifts 
and with different optical activity classifications. Neu- 
tral and ionized gas tracers are affected, as are lines with 
different rest-frame wavelengths, intrinsic optical depths 
and critical densities. This eliminates most of the pro- 
posed explanations for the [C II] deficit, leaving the de- 
pendence on the Lfir/Myi^ ratio as the most plausible 
dominant factor driving the deficits. We have assumed 
in this paper that ipiR is mostly due to star formation. 
This might not be true for some galaxies like Mrk 231, 
but even an AGN contribution of 50% to ipiR will fall 
short of explaining the factor of 10 decrease in some of 
the line to FIR ratios. 

In the next section we show with quantitative modeling 
that high ionization parameters can indeed account for 
most of the observed line deficits. The key factor is that 
in dusty star forming regions the fraction of luminosity 
absorbed by dust particles in ionized as well as neutral 
gas increases with U. 

4. COMBINED HII AND PDR LINE MODELS 



lAbel et ahl (|2009|) used t he spectral synthesis code 
Cloudv (jFerland et aLlfTOOl ) to model the [C II] /FIR ra- 
tio in clouds exposed to different radiation fields. They 
found that they could reproduce the [C II] line deficits 
observed in some local LIRGs and ULIRGs, as well as 
their warmer IRAS S'eo/im/'S'ioo/im colors, by increas- 
ing the value of the ionization parameter at the surface of 
the clouds. As U increases, the HII region is extended to 
higher Ay into the cloud. As a result, a larger fraction of 
the UV photons are absorbed by the dust in the ionized 



4 



Gracia-Carpio et al. 




Figure 3. Cloudy model predictions for the observed line to FIR continuum ratios as a function of the hydrogen density and the ionization 
parameter at the illuminated face of the cloud. Red (grey) areas indicate the approximate range of values in the line to FIR continuum 
ratios observed in galaxies with Lfir/Mhj > 80 Lq Mq~^ (< 80 Lq Mq~^). The dashed contour indicates the values of the ionization 
parameter and hydrogen density that can simultaneously explain all the line to FIR ratios in galaxies with Lpir/Mh, < SO Lq Mq~^ . 
The arrows show the changes to the line to continuum ratios when the ionization parameter is increased by an order of magnitude. 



region and reemitted in the form of infrared emission. 
The net effect is that the fraction of UV photons available 
to ionize and excite the gas is reduced at high U, decreas- 
ing the relative intensity of the fine st ructure lin es com- 
pared to the FIR continuum (see also lWtllT99l) . bmce 
the number of UV photons per dust particle increases 
with U, the dust temperature (and the S'eo/im/'S'ioo/im 
color) increases. 

We have extended the analysis in lAbel et all ()2009[ ) to 
include all the fine structure lines observed with PACS. 
We used Cloudy to simultaneously calculate the emission 
from the molecular, neutral and ionized gas assuming 
pressure equilibrium in the clouds. The ionizing source 
was re presented by a non-LTE CoStar stellar atmo- 
sphere (jSchaerer fc de Koted[T99l with = 36000 K, 
and the calculations were stopped at Ay = 100 mag 
{Nh = 2 X 10^3 cm-2 for Av/Nh = 5 x 10^22 magcm^). 
To facilitate th e comparison with previous studies 
(|Kaufman et all Il999; Abel et al. 2009), we used the 
.as-ph ase elemental abundances from lSavage fc SembachI 
19961) : C/H = 1.4 x 10"'', 0/H = 3.2 x 10"^, N/H 
= 7.9 X 10"^. In our model the spectral energy distribu- 
tion of the ionization source is fixed. Thus, for a given 
density, an increase of the ionization parameter results 
in a proportional increase in the number of UV photons 
at all energies. 

In Figure [3] we show the model predictions for the line 
ratios represented in Figure [21 Models with U ~ 10"^- 
10~^'^ and hydrogen densities ~10^'^-10^'^ cm"^ can si- 
multaneously explain the average line to FIR ratios mea- 
sured in galaxies with low Lfir/Mhj values. The ob- 
served [C II] , [N II] and [O I] line deficits in galaxies with 
Lfir/Mh^ > 80 Lq Mq~^ can be easily understood if 
their ionization parameters are at least an order of mag- 



nitude higher, consistent with the idea that the ionization 
parameter is approximately proportional to Lfir/A/h2- 
However, we cannot explain with this particular model 
the observed line deficits in the [O III] and [N III] high 
ionization lines. These lines are produced closer to the 
ionizing source, where the effect of the dust opacity on 
the formation of the lines is less important. This may 
indicate that the ionization parameter is not the only 
property that varies as L-piji/ Mn^ increases. In future 
papers we will explore the full parameter space in these 
models, varying the ionization source (SED), the stop- 
ping criteria and the metallicity of the gas. We will also 
study the possibility of multiple components contribut- 
ing to the global far-infrared spectrum of these galaxies, 
including partial coverin g of the high ionization regions 
(e.g.. . Fischer et al.|[20Tol ). 

5. DISCUSSION AND CONCLUSIONS 

As discussed in the introduction, galaxies can expe- 
rience short periods of intense star formation, normally 
associated with the merging of massive galaxies, in which 
the molecular gas is transformed into stars with high ef- 
ficiency. A good way to investigate this enhancement 
of the star formation efficiency is to study the relation 
between the molecular gas surface density, Sjja , and the 
star formation rate surface density, Ssfr- Recent studies 
of the global infrared continuum and molecular gas prop- 
erties of galaxies in the local and high redshift Universe 
suggest that, independent of redshift, luminous mergers 
and non- or weakly-interacting star-forming galaxies fol- 
low two separate Kennicutt-Schmidt relations (Ssfr oc 
Sj^^) with similar exponents N ~ 1.1-1.2, but different 
normalizations ijGenzel et all 120101: iDaddi et al]l2010D . 
The explanation for this offset in the star formation rela- 
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tions is still not clear, but it may indicate that, in addi- 
tion to Sh2: other galaxy properties, like the dynamical 
time scale, may play a significant role regulating the star 
formation efficiency of the gas. 

The Lpin/AfHa value that separates galaxies into these 
two types of star formation {^100 Lq Mq~^) is remark- 
ably similar to the point where we start to find a decline 
in the fine structure line to FIR ratios. This indicates 
that the average properties of the neutral and ionized 
gas are significantly different in these galaxies, broadly 
consistent with the scenario of a highly compressed more 
efficient star formation, creating largely enhanced ioniza- 
tion parameters that manifest themselves in lower line to 
continuum ratios. Our detection of a universal drop of 
the line to continuum ratio at about the same Lfir,/-^^H2 
value where Genzel ct al. (2010) and Daddi ct al. (2010) 
claim a transition to a more efficient star formation mode 
gives in turn credence to this interpretation. 

In the next two years Herschel will increase the num- 
ber of galaxies with fine structure line observations in the 
far-infrared, especially in objects with high Lfir/-^^H2- 
The detailed study of the line ratios, combined with the 
use of spectral synthesis codes like Cloudy, will provide 
valuable information about the ionizing sources, possible 
AGN contribution to Lfir, far- infrared optical depth, 
and the physical properties and metallicity of the gas. 
This information will be critical to understand the driv- 
ing mechanisms involved in these modes of star forma- 
tion. 
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(Italy); lAC (Spain). This development has been sup- 
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